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Current auto- and cross-correlations are studied in a system of two capacitively coupled quantum 
dots. We are interested in a role of Coulomb interaction in dynamical correlations, which occur 
outside the Coulomb blockade region (for high bias) . After decomposition of the current correlation 
functions into contributions between individual tunneling events, we can show which of them are 
relevant and lead to sub-/supper-Poissonian shot noise and negative/positive cross-correlations. The 
results are differentiated for a weak and strong inter-dot coupling. Interesting results are for the 
strong coupling case when electron transfer in one of the channel is strongly correlated with charge 
drag in the second channel. We show that cross-correlations are non-monotonic functions of bias 
voltage and they are in general negative (except some cases with asymmetric tunnel resistances). 
This is effect of local potential fluctuations correlated by Coulomb interaction, which mimics the 
Pauli exclusion principle. 

PACS numbers: 73.23.-b,72.70.+m,73.23.Hk,73.63.Kv 

I. INTRODUCTION 

Recently, McClure et alP and Zhang et alP performed measurements of current noise auto- and cross-correlation 
in a system of two capacitively coupled quantum dots (2QD). It is first measurement showing some aspects of current 
correlations in 2QD, as anti-bunching and bunching of scattered electrons. The 2QD system is simplest, in which one 
can see a role of Coulombic interactions on current-current correlations. 

In a system of noninteracting electrons the Pauli principle is essential in scattering process and it can lead to anti- 
bunching P This effect in a multiterminal geometry, such as in the Hanbury Brown and Twiss experiment, ^ can be 
manifested in negative cross-cor relatio n between scattered electrons. The phenomena was studied theoretical!^ 6 -^ and 
observed in several experiments! 8 * 9 * 1 ^ In nanostructures anti-bunching is responsible for reduction of auto-correlation 
shot noise Si j below the Poissonian value Sp = 2el, where / is the average current and e is the charge of an 
electron! 11 * 12 * 13 * 1 ^ 5 -! For sequential tunneling through a single quantum dot (QD) a maximal reduction of the shot 
noise can be to the value 1/2 x 2e/! 15 ' 16 1 17 1 1 ^ 

An enhancement of shot noise above the Poissonian value was observed in a resonant tunneling diode by Iannaccone 
et alP^and in a quantum well by Kuznetsov et al. 20 in the region of negative differential resistance. The enhancement 
of shot noise was interpreted as a result of change in the density of states, which causes more state to be available for 
tunneling from cathode. Strong back-scattering, at cu rrent pinch-off (close to the Coulomb blockade border), can be 
also responsible for an enhancement of shot noise l 2 ^ 2 ^ In the Cou lomb blockade region inelastic cotunneling processes 
can affect transport and lead to super-Poissonian shot noise 1 \22]M 

Recently, Wu and Yipf^ and Rychkov and BiittikeiJ^ applied a phenomenological Langevin formalism in circuit 
modeling to calculations of current cross-correlations in a 3-terminal system. An interesting result is that current 
cross-correlations are always positive in a macroscopic classical 3-terminal with a fluctuating current in an input 
electrode. It is in contrast to a microscopic picture of scattered electrons, which in general shows anti-bunching with 
negative cross-correlations ! 5 * 6 ^ However, there is an exception when inelastic scattering occurs. Texier and Biittikei^ 
considered a multiterminal system with inelastic scattering in an additional electrode, in which the current was kept 
equal to zero and which caused voltage fluctuations. They showed that inelastic scattering can lead to positive cross- 
correlations, whereas for quasielastic scattering correlations remain always negative. Just recently Oberholzer at 
al! 2 ^ confirmed experimentally these predictions in a device with fractional quantum Hall effect edge states, in which 
interactions between current carrying states and a fluctuating voltage were controlled by an external gate voltage. If 
scattered electrons are entangled, the situation can be different, one can see bunching or anti-bunching. Burkard et 
al! 2 ^ showed that for singlets the auto-correlation function is twice as large as for independent electrons and is reduced 
to zero for triplets. 

Cur rent corre lations in the system of two capacitively coupled quantum dots were studied theoretically in several 
papersP ^ 30 1 31 1 32 1 3 ^ Our studies on a ferromagnetic single electron transistor (fSET)PMIl are closely related, because 
there are two channels for electrons with opposite spins a =j, j and Coulomb interactions are included as well. For 
the first time we pointed on dynamical Coulomb blockade effect between the both channels, which is responsible for a 
reduction of current and leads to the super-Poissonian shot noise. Positive cross-correlations were predicted by Cottet 
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et al! 35 * 36 * 3 ^ on a similar system with three terminals and 3QDs. They also analyzed time evolution of tunneling 
events and presented bunching. Dynamical Coulomb blockade and bunching were studied by Gustavsson et al! 38 | S9 1 
in time-resolved measurements of electron transport through a multilevel QD. 

In this paper we want to study a microscopic nature of dynamical current corre lations, especially a role of Coulomb 
interactions, in the 2QD system. Some problems were mentioned earlier™^ for example, a role of charge and 
polarization fluctuations and their contribution to the current correlation function. We expect that for a symmetric 
case these contributions may compensate each other, but they can be pronounced in an asymmetric case. Two 
different situations will be studied, with a weak and strong inter-dot coupling. We will show that both the cases 
have different charging diagrams in the space of biases applied to the both channels. For the strong coupling electron 
transfer in one of the channel should be strongly correlated with charge drag in the second channel. In order to 
see correlations between individual tunneling events we decompose correlation functions and analyze contributions 
leading to the sub-/super-Poissonian noise as well to positive cross-correlations. The experiments^] were performed 
for small biases and near degeneracy lines separating different Coulomb blockade regions. The present studies are 
more general. We analyze the Schottky contribution, which can be large close to the Coulomb blockade border due 
to backward scattering. However, our main interest focus on dynamical correlations for electronic transport outside 
the Coulomb blockade region when the bias window is large. 

The paper is organized as follows. In the next section II we describe the 2QD model, derive general formulas 
for charging energy, current and correlation functions. The section III presents charging diagrams and allowed 
tunneling processes in the space of bias voltages. The charging diagrams are background for further analysis of 
correlations between individual tunneling events presented in the section IV. We show there differences of auto- and 
cross-correlations for weak and strong inter-dot coupling. The correlation functions are non-monotonous with bias 
voltages, which means that there are several competing correlation processes. Our decomposition approach shows 
which contributions are relevant. For example, current cross-correlations are in general negative (except special 
situations with asymmetric resistances), but this anti-bu nching effect is due to potential fluctuations (it is not due 
to the Pauli exclusion principle as in many nanostructures^ 12 * 13 * 14 * 1 ^ ) . Main results of the paper are summarized in 
the final section V. 



II. DESCRIPTION OF THE MODEL AND GENERAL DERIVATION OF TRANSPORT QUANTITIES 



A. Model 



We consider a system which consists of two large capacitively coupled quantum dots (top and bottom: £ = t, b) with 
a continuous electronic density of states. Each quantum dot (QD) is connected with the left (a = L) and the right 
(a = R) electrodes through the tunnel junctions which are characterized by the tunnel resistances Ri a and capacitances 
Cia (see Fig. [I]). The strength of the inter-dot Coulomb interaction is modeled by the capacitance Ci n t. To this 4- 
terminal system, the bias voltages V t R, VbR are applied to the right electrodes, while both the left electrodes are 
grounded (asymmetric bias). Usually additional gates are applied to the QDs, which can change their electrochemical 
potentials and number of electrons. In order to simplify the analysis these gates can be omitted. The corresponding 
resistances of the tunnel junctions are assumed to be much larger than quantum resistance Rq — h/2e 2 so the 
electronic transport is dominated by incoherent, sequential tunneling processes™ Since our research is focused on 
dynamical processes in high voltage regime (outside the Coulomb blockade), higher order processes (e.g., cotunneling) 
are neglected. Cotunneling can be relevant in the Coulomb blockade region, for tunneling currents as well as super- 
Poissonian shot noise) 3 -^ 31 * 32 * 33 * It is also assumed that the temperature T is low which means that charging energies 
of the individual dots Ef 1 and an energy of the inter-dot Coulomb interaction E£% t are much larger than a thermal 
energy k B T, i.e., Ef, E$ t > k B T. 



B. Charging energy 



The voltage drops on the tunnel junctions as well as between both dots (i.e., on the capacitance Ci nt ) can be 
calculated from Kirchhoff 's laws: 



V tL = ~^r n t + —7- nb + c2 V t R H ^ — VbR (I) 



(2) 
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V b d [ op = ^fn t + ^fn b + ^^V tR + ^V bR (3) 
V b t P = V bR -V b d [ op (4) 

T rdrop ^(CbL + CbR) e(C t L + C t R) , (CbL + CbR)C t R j j- (C t L + C t R)C b R , . 

V int = ^2 U t C~2 n *> + ^ V tR C~2 ( 5 ) 

where an electron charge e < 0, C t = C tL + C tR + C int , C 6 = C 6L + + C int , C 2 = C t C h - C 2 nt . The charging 
energies Ef = e 2 C b /C 2 , E£ h = e 2 C t /C 2 and E£% t = e 2 C int /C 2 . The inter-dot energy E^ t describes the change in 
the energy of one dot when an electron is added to the other dot. 

Transfer rate Ff a through the junction £a can be found from a changes of the free energy of the whole system when 
an electron is added (s = +) (or extracted: s = — ) to (from) one of the quantum dots. The free energy consists of the 
potential energies of the electrodes and the electrostatic energies of the charged capacitors and for the system under 
consideration has a form: 

F[n u n b , n t Li n t R, n bL , n bR ) = -QtReiVtR - QbReiV b R 

c tL (v t d rf + c tR {v t t p ? + c bL (v b d D 2 + c bR {v b d R p ? + c mt (vf;n 2 ] /* , ® 

where QtRei = QtR + en t R and Q b R e i = Q bR + en b R are charges on the top and bottom right electrodes, respectively. 
ntR (nbR) is a number of electrons transferred from the right electrode to the top (bottom) QD and charges on the 
capacitors Ce a are Q £a = Ci a Vf™ p (£ = t, b; a = L, R). The changes of the free energy due to electron tunneling 
through the junction £a can be written as a differences between the free energies of the initial Fi and the final Ff 
states 

AF/ L (n t , n b ) = Fi(n u n h , n tL , n bL , n tR , n bR ) - F f (n t + si, n h , n tL + si, n bL ,n tR , n bR ) , (7) 
AF^ L (n u n b ) = F i (n u n b ,n tL ,n bL ,n t R,n bR ) - F f (n u n b + sl,n tL ,n hL + sl,n tR ,n bR ) , (8) 
AF£ R (n u n b ) = Fi(n u n h , n tL , n bL , n tR , n bR ) - F f (n t + si, n b , n tL , n bL , n tR + si, n bR ) , (9) 
AFjj R (n t ,n b ) = Fi(n t , n b , n tL , n bL , n tR , n bR ) - F f {n u n b + si, n tL , n bL ,n tR , n bR + si) , (10) 

where s = ± for an added or extracted electron. An electron can be transferred through the junction lot when the 
corresponding free energy difference AFg a (n t ,n b ) > 0. 



C. Current and correlation functions 



One can calculate the tunneling rates through the junctions using the method developed for a single QD with a 
continuous electronic density of states (DOS) (see, for example, 40 ). For our system, one can obtain the following 
tunneling rate for transfer of an electron through the junction £a: 

V s , \ \ AFj a (n u n h ) 

l £a [n u n b) - e2R ^ x _ exp[ _ AFL{nunb y kBT] > W 

where the tunneling resistance R £a = h/ (Ane 2 \M £a \ 2 D £ D a ) and D £ is DOS in QD, while D a is DOS in the electrode. 
It is assumed that the transfer matrix element M £a as well as DOS are constant around the Fermi energy, and the 
resistances R £a are parameters of the model. 

In the stationary state the average currents flowing through the junction £a can be found from the formula 

It a = (S La -S Ra )(I+-IZ x ), (12) 

where 8 La, $Ra are the Kronecker's deltas, If a = e ^Z nt Uh Tf a (n t , n b ) p(n tl n b ) is the current flowing to/from the 
QD (s = ±), and the probability p(n t: n b ) describes the system in the steady state which contains n t and n b excess 
electrons on the top and the bottom QD, respectively. The probability p(n tl n b ) can be found from the master 
equation: 

dv ^ nt ^ b)t ^ = ^ [F s ta (n t - si, n b )p(n t - si, n 6 ; t) + rg a (n £ , n b - sl)p{n t , n b - si; t)] 

a,s 

-p(n u n 6 ; t) ^ [r? a (n t , n b ) + T 8 hOL (n u n b )} (13) 
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with the left hand side equal to zero. 

To analyze fluctuations in the system we extend the generation-recombination approach^ for multi-electron channels 
by a generalization of the method developed for spinless electrons in a SET.^ In calculation of the current-current 
correlation functions, we also include the self-correlation terms as wellP^ According to this procedure the auto- and 
cross- current- current correlation functions are 

Se a ,e'a'(w) = $w5oLa'Sfu h + S£ a ,£'a'( u )i (14) 
where the Schottky term (the frequency independent part for uo — > oo) is given by 

S$' h = 2e(I+ a + /,"), (15) 

and the frequency dependent part 

S? ai/ , a ,M = 2e 2 (S La - 5 Ra ){5 La , - 5 Ra ,) £ D<<^{ui) . (16) 

n' t ,n' h \n t ,n h 

Here, we denoted 

Df a % tn \ui) = J2 SS ' { r «K-<) G n , <tntnb {uj)[5 u Tl a {n t - sl,n b )p(n t - sl,n b ) 

s,s' 

+ s beT ^(nt, n b - si) p(n u n b - si) ] 
+^ a ( n ^ n 6)^n;n;,n t n b (-^)[^r| w (n^ - s'l,n b )p(n t - s f l,n h ) 

+8wTi, a , (rit, n b - si) p{n u n b - s'l) } } (17) 

as contributions to the dynamical part of the correlation function for various tunneling events in the space of the charge 
states (n t ,n b ). Later we will analyze the components D™^ n i b /™ / tnb , in order to show correlations between tunneling 
processes. The elements of the Green function G n ' tn ^ ntnh {uo) = {iuol — M)~} n , ntfib and the matrix M is constructed 



from the right hand side of the master equation (13). Here, we use the symmetrized formulas [Eqs. (|16|)-( 17)] for the 
current shot noise. Recently, an asymmetric formula was used to studies frequency-dependent asymmetric features.^ 
Our analysis is confined to the zero-frequency limit, and therefore, the symmetrized approach is justified. 
We can write also other correlation function for any quantity X and Y in the space of states (n tl n b ) 

Sxy{u) =4 ^2 x ( n ^K) G n f t nln t n b (^)y(n tl n b )p(n t ,n b ) . (18) 

n' t ,n' h \n u n h 

For example, in further studies we will calculate the volt age- volt age correlation function Sy t y t , where V t = (E^ h n t + 
E^ t n b )/e and it is a potential difference at the top QD when an additional charge is added to the system. The 
function Sy t y t describes potential and charge fluctuations during current flow. 



III. CHARGING DIAGRAMS 



Before solving the master equation (13), we should first to determine [from Eq. ( pTj )] the tunneling rates T s £a (n t , n b ) 
in the two-dimensional (2D) space of charge states {n t ,n b } for given bias voltages V t R and V b R. If AF^ a (n tl n b ) < 
then Tf a (n t ,n b ) ~ exp [—\AFg a (n t ,n b )\/kBT] /(e 2 R£ a ) and the tunneling process is exponentially suppressed in the 
low temperature limit (T — > 0). Otherwise, when AF^ a (n t , n b ) > 0, the tunneling process can be relevant for transport. 
In practice one can restrict considerations to the {n t ,n b } space around the ground state (0,0) with a small number 
of available states and tunneling processes relevant for transport. Of course, the space of available states depends on 
the bias voltage Vi a and increases with an increase of Vi a . If one wants to consider temperature effects, the space 
{n tl n b } should be enlarged and neighboring charge states with thermally activated processes should be taken into 
account. 

Here, we perform an analysis of the available states and relevant tunneling processes at T = for a weak and 
strong inter-dot coupling case. First, the weak coupling case is considered, for which one can expect that transport is 
similar to the case of two independent QDs and one can easily understand modifications introduced by the inter-dot 
coupling. We consider the system with symmetrical couplings to the electrodes C t L = C t R = C b L = C b R, in which the 
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charging energy of the both QDs is Ef 1 = E^ h . This case is simpler to presentation. It is no problem to generalize 
the analysis for different couplings da- The available charge states {nt,n&} and tunneling processes participating 
in currents are presented in Fig. [2] in the (V^V^r) space. The diagram was obtained for dnt = Ci a = 1 aF. The 
rhombic area in the middle of the Fig. p[a) corresponds to the Coulomb blockade (CB) region, where all TJ a (0,0) are 
suppressed and electrons can not leave the state (0,0). With increasing voltage V t R (or V^r) some tunneling processes 
become allowed, because the corresponding free energy differences become positive, AF/ a (0,0) > 0. New charge 
states (n t ,rib) become available and the currents begin to flow. Threshold voltages can be determined from the 
equation AF| Q ,(nt, rib) — and they are presented in Fig. [2^a) as (green) solid lines. The threshold lines divide the 
(VtR, VbR) space into regions {a, . . . , n} with particular allowed tunneling processes {Tf a (n t , rib)} presented below in 
Fig. [2^b). Let us analyze as an example transport for the case of small VbR ~ 0. With an increase of V t R one can 
reach above a threshold voltage V th the region 6, where r^(0, 0) becomes relevant and electron can leave the top QD 
through the right tunnel junction. In this case only two charge states (0,0) and (—1,0) participate in transport. In 
this region one can expect similar current and shot noise characteristics as those for independent QD. The tunneling 
current is flowing only through the top QD while the bottom QD is still in the Coulomb blockade regime. On the 
other hand, in the region / only two charge states (0, 0) and (0, 1) participate in transport so the tunneling current 
is flowing through the bottom QD while the top QD is in the CB regime. In the region a the voltage VbR is higher, 
r^ L ( — 1,0) becomes relevant (tunneling through the bottom left junction is allowed) and four charge states (0,0), 
(—1,0), (—1,1), (0, 1) participate in transport. On the left hand side of Fig.|2|a) we have VtR < VbR and the situation 
is symmetric to that one on the right hand side. For higher voltages, outside the regions {a, . . . , n}, additional states 
(nt,rib) are available for the tunneling events. 

The charging diagram and the relevant tunneling processes for the strong inter-coupling case (Ci nt — 3 aF) are 
presented in Fig. [3j We have found that the ground state is always (0, 0) in the CB region, with a global minimum 
of the system. However, a set of local minima appears at (n, — n) (with n = ±1,±2, ...). These polarized states 
(n, —n) can be only reached in CB through thermally activated neighboring states (n, — n ±1), (nil, — n). Because 
the corresponding free energy increases with the polarization P = n t — rib, thus these states are less stable for large 
P. The system evolves between the local minima in order to reach the ground state, and transition times are shorter 
between the metastable states with a large polarization P. One can expect that the polarized states and fluctuations 
of polarization play a relevant role in currents and shot noise. For the case dnt = 3 aF the local minima appear 
for the polarized states (1,-1) and (—1,1), and all transfer rates r| a (l, —1) and r| a (— 1, 1) are exponentially small 
in the CB region. Comparing the diagrams in Fig. [2Va) and|3ja) one can see that charged regions shrink with Ci nt . 
The CB region is smaller as well as charged regions {a, ...,n} above V th . For example, in the region b (at VtR > V th 
and VbR « 0) all four states (0,0), (—1,0), ( — 1, 1), (0, 1) participate in transport for dnt = 3 aF, whereas in Fig. [2] 
only two states (0,0) and (—1,0) are available. For the strong inter-dot coupling the bias voltage V t R opens only the 
transfer (0,0) — > ( — 1,0), since all other tunneling processes are allowed. One see a close loop of electron transfers 
between the charge states (0,0) -> (-1,0) -> (-1, 1) -> (0, 1) -> (0,0). 

IV. ANALYSIS OF CURRENTS AND SHOT NOISE 
A. Weak inter-dot coupling 

In this section we want to study dynamical current-current correlations. Therefore, we analyze contributions of 
different tunneling processes Tf a (n t , rib) m the (n tl rib) charge state space to currents Ii a and shot noise Si a ^ a f. Their 
voltage characteristics depend on the charging energies, on AFj> a (n t , n^), as well as on the tunnel resistances Re a - Fig.H] 
presents results of numerical calculations for the weak inter-dot coupling Ci nt = 1 aF. We plot maps in the (VtR, VbR) 
space for the current I t L through the top QD, the Fano factor F t L = StL,tL(w = 0)/2eI t L (for the auto-correlation 
current-current function in the top QD at the frequency uo = 0), the cross-correlation functions StL,bL — $£L,6l(^ = 0) 
between the currents in the top and the bottom QD, and the Fano factor FbL — SbL,bh(u — Q)/2eIbL at the bottom 
QD. The current IbL is not shown, because its characteristic is very similar to I t L (one can rotate its map by 
90 degrees). The tunnel resistances are assumed to be symmetric for the top QD (RtL — RtR — 1 M.Q) and 
for the bottom QD (RbL — RbR — 50 MQ). We have analyzed the characteristics for other resistances (also for 
RtL = RtR = RbL = RbR = 1 M^), but the super-Poissonian shot noise is more pronounced in the case presented. 
Because temperature only smears the characteristics close to the threshold voltages V th , so it is irrelevant for our 
studies, and we take T = in the calculations. We do not analyze the CB region, where other (cotunneling) processes 
can be relevant, and where they can dominate over contributions of sequential tunneling to currents. 

Let us analyze the Fano factor F t L for the top QD, which is presented in Fig. ^b). In the region b (for V t R > V th 
and VbR w 0) transport through the top QD is not disturbed by the bottom QD, and the shot noise is in the sub- 
Poissonian regime (with F t L < 1) as expected. For larger positive values of VbR (in the region a, c, or d) the Fano 
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factor F t L > 1 and decreases below the unity with an increase of V t R. In this case the charge states (0,0), (—1,0), 
(—1,1) and (0, 1) participate in transport. Charge fluc tuations in the bottom QD can be large and can lead to the 
dynamical Coulomb blockade (DCB) effect \MMMMM1 The effect manifest itself in an enhancement of the shot noise 
to the super-Poissonian regime (with F t L > 1). It is seen in the frequency dependent part of the auto-correlation 



function, S^ LtL {uo) [see Eq. (16)], which at uo = is a sum of all contributions of current fluctuations with various 
relaxation times p3M This function is then positive, since many tunneling processes show bunching of electrons and 
their contributions are positive. Fig. [4jb) shows that the super-Poissonian shot noise occurs also in the regions e and 
9- 

In order to understand a role of various tunneling processes in the shot noise we analyze all components 
D^l n t b £ ntnb , which contribute to the dynamical part of the function S% L tL [see Eqs. il6Ml7l]. The results are presented 
in Fig. |5^a), which is a cross-section through the regions a, 6, c and d in Fig. [4pD). The auto-correlation function 
D~ L describes correlations of tunneling processes between the charge states ( — 1,0) — > (0,0). In the region b the 
function is negative, and therefore, the shot noise is sub-Poissonian. But the function increases and becomes positive 
in the regions a and c. A more pronounced increase is seen for the component D^tL 1 • ^ n contrast, the cross- 
correlation function D^®^ 11 is negative in all the regions. It means that the tunneling processes (—1,0) — > (0,0) 
and (—1, 1) — > (0, 1) are anti-bunched. This contribution, however, does not compensate D^®^ 10 and D^^ 11 , and 
therefore, FtL > 1 in the regions a and c. Our previous studies^^ showed that the super-Poissonian shot noise is 
due to activation of charge fluctuations in the system. For the present model we have the same situation. Fig. |5|b) 
presents the plots of voltage correlation functions Sv t ,v t and Sy b ,v b in the top and the bottom QD. It is seen a large 
enhancement of voltage and charge fluctuations in the regions a, c and d in the bottom QD (see the blue dash curve). 
Moreover, these plots are very similar to the auto-correlation components D^l n t b £ ntnb , even small kinks at V^r = —36 
mV are similar in both the figures. Voltage fluctuations in the bottom QD are very strong, two orders of magnitude 
larger than those in the top QD. 

The cross-correlation function StL,bL is presented in Fig.[4^c). We do not normalize the cross-correlation functions, 
because the function is a sum of two contributions proportional to some components of I t L a nd I^l, and therefore, 
there is not a good normalization factoJ^ (neither I tL , nor I bL , nor y/TtiJbL can be use d 35144 1 ). One see in Fig.gc) 
that the function S t L,bL is negative in the right-upper and the left-lower quarter of the map, when the currents I t L and 
IbL have the s ame direction. T his dependence is typical for cross-correlation in scattering of fermions, which exhibit 
anti-bunchingPI 1 * 12 * 13 * 14 * 15 * 43 ! However, in our situation electronic transfer is separated in the top and in the bottom 
part of the device, so anti-bunching is not fermionic nature. We show later that correlation of both the electronic 
channels is due to Coulombic interaction and potential fluctuations. 

In contrast to FtL, the current noise F^l in the bottom QD, is generally the sub-Poissonian type, except with a 
narrow stripe around V\>r = for high V t R [see Fig. [4^d)] . The difference between the maps in Fig. [4jb) and|4^d) 
results from the DCB effect. An electron traversing the bottom channel spends a long time at the bottom QD and 
the charge fluctuations at the top QD are relatively fast. Therefore, transport in the bottom channels depends on an 
average potential in the top QD and the shot noise F^l is very little sensitive to individual tunneling events in the 
top channel of our device. 



B. Strong inter-dot coupling 

Maps in Fig. [6] show transport characteristics ItL, FtL, StL,bL and F^l in the bias space (VtR, VbR) for the strong 
inter-dot coupling. In this case the maps are different than for the weak coupling in Fig. [4j because the charging 
diagram is different (compare Figs. [3] and In the regions {a, all four charge states (0,0), (—1,0), ( — 1,1), 

(0,1), contribute to transport; whereas the states (0,0), (1,0), (1,-1), (0,-1) are active in the regions {h : 
Tunneling processes and their contributions to transport are different than in the case of the weak coupling. Let us 
analyze the Fano factor F t L presented in Fig. (6^b). In the region b (for V tR > V th and V bR « 0) the shot noise is 
super-Poissonian, but one see also a dark area (with F t L < 1) inside this region. It means that there are various 
competing tunneling processes and their role changes with the bias voltage. An interesting situation is in the region 
/ (for VbR > V th and small Vt R ). In the middle F t L shows a peak, while in Fig. Qb) the shot noise decreased below 
the unity. One see the dark areas in Fig. [6^b) close to the borders with the regions e and g. 

We have analyzed several cross-sections of the maps in Fig. [6j One of them is in Fig. [7J which presents F t L in the 
region b at V bR = 0. We also plotted the dynamical part of the shot noise, F£ L = S^ L tL /2eI t L, and its components 

Dtl n t b £ ntnb /2eI t L- For small voltages V bR the auto-correlation components D~^^ 10 and D^ 1 ^ 11 dominate over 
the cross-correlation one DtL^tL 11 ^ so ^tL tL > and the corresponding F t L is super-Poissonian. It is seen that 
the function D~^^ n always decreases with V tR . The auto-correlation components have different dependencies, 
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and these two components compete each other. The function D tL 1 ^ L 10 increases, whereas D tL 1 ] , L 11 decreases. For 
small bias voltage tunneling processes (—1,1) — > (0,1) dominates in the super-Poissonian shot noise, while auto- 
correlations (—1,0) — > (0,0) become dominating for large V t R. In the intermediate voltage range F t L < 1, because 
the cross-correlation component plays a dominating role. 

In the Fig. [8] we present F t L and its relevant components at VbR = — 50 mV for a cross-section through the regions 
/, e and d. The strong enhancement of the F t L in the region / at small VtR is due to the Schottky term Sf£ h . This 
situation is in contrast to that one for VbR = presented in Fig. [7| where S^£ h /2eI t L = 1 in whole region b (with 
exception of a small region close to the threshold voltages, where strong thermal current fluctuations occur). From 
the plots in Fig. [8] one see that the contribution Sf£ h to the total Fano factor decreases with increasing voltage V t R, 
while a role of the dynamical part S% L grows. This leads to the reduction of the Fano factor below Poissonian value for 
voltages 14.5 mV < V t R < 24 mV (close to the border between the region / and e). For higher voltages S% L tL become 
positive and the Fano factor is super-Poissonian. One see from Fig. [8] that an enhancement of the F^ in the region e 
is caused by auto-correlation processes F>^ L 1 ^ L ~ 10 (red dot curve). The border between regions / and e is a transition 
line between the sub- and the super-Poissonian shot noise, where the tunneling process T^ R (— 1,0) disappears while 
the process r^(0,0) becomes activated. 

The cross-correlation S t L,bR [see Fig. |6jc)] is similar to that one in the weak coupling case, which shows anti- 
bunching. However, this anti-bunching process is not fermionic origin, exchange of electrons is impossible between 
the both QDs. In our system strong Coulombic repulsion occurs between charges transferred through the top and 
the bottom part of the device. This process mimics the Pauli exclusion principle, and therefore, electronic transfers 
through the top and the bottom QD should be anti-correlated and the function S t L,bR < 0. We make a cross-section 
through regions a, b and c at V t R = 50 mV for the current IbL and S t L,bR [see Fig- [9|a)]. We expected that the 
cross-correlation function StL,bR is proportional to the current —IbL, at least in the region b for small bias. However, 
the plot is different - its shows a plateau around VbR ~ and a sharp increase, which is proportional to —IbL- It is 
clear that the current IbL and the cross-correlation function S t L,bR may change their dependencies in the regions a 
and c, because new tunneling processes are activated. In order to understand the dependence of S t L,bR i n the region 

b we plot its components D™£ l £j™ tnh in Fig. j^b). Notice that at VbR = all the components are non-zero, but they 
compensate each other and the sum is S t L,bR = 0. It is clear that D~ L 1<d hL ~ 10 + D^^f 1 = at VbR = 0, because 
these components describe correlations between the tunneling event (—1,0) — > (0,0) through the left top junction 
and the forward and the backward tunneling event through the left bottom junction (between the charge states 
( — 1,0) — > (—1, 1) and (0, 1) — > (0,0), respectively). Similarly, one has for D~ L l \ L ~ 10 + D~ L ll h ^ 1 = 0. However, one can 

see that at VhR = —3.55 mV the components D^l^l 10 = and D^l^l 1 — 0, while D~^^ 10 = and D^^f 1 = at 
VbR = —3.73 mV. This compensation is due to advanced and retarded correlations, which occur for any component 

^th'bL^ I see |^| an d terms corresponding to u and — u}. 42 Both the contributions, advanced and retarded, 
compensate each other either at VbR = —3.55 mV or -3.73 mV, respectively. This effect is also responsible for the 
plateau seen in S t L,bR- Outside the plateau region the dominating components are D^®'® 1 and D~ L X \^ 10 (for negative 
and positive VbR, respectively) - see the pink long dash and the blue short dash curves in the bottom part 
Let us come back to the maps in Fig. [6] and analyze the Fano factor FbL in the bottom QD shown in 
It is the sub-Poissonian type in the most regions, but the function has an interesting feature in the region b and z, 
where FbL > 1- The shot noise is due to charge fluctuations in the top QD. When an electron leaves the top QD, 
(0, 0) — > (—1, 0), the system becomes unstable, therefore, an electron is attracted to the bottom QD, ( — 1, 0) — > ( — 1, 1), 
through the left or the right tunnel junction. The system is in the metastable state (—1,1). The next tunneling event 
in the top QD, (—1,1) — > (0,1), leads to pushing out the electron from the bottom QD, (0,1) — > (0,0). This is a 
pumping effect - a charge injected into (or ejected from) the top QD leads to ejection (injection) of a charge at the 
bottom QD. One can expect that the shot noise SbL,bL should be proportional to the tunneling current 7 t L- It is easily 




to prove it, calculating these quantities in the region b for VbR = 0. The results are presented in Fig. 10 The main 



contribution to the shot noise [Eq. 14 is due to the Schottky term [Eq. 15 , which is now S^£ h = 4e7^. Since for 



VbR = the current flowing into and from the bottom QD are equal iX = L, the total current IbL = l£r — I bL = 0. 



Therefore, Fig. 10 presents the auto-correlation function SbL,bL, instead the Fano factor FbL = SbL,bL/2eIbL- In this 

mula 

T+ _ r t - fl (o,o)r+(-i,i)r+(-i,o)r fc -(o,i) 

(20) 



case one can obtain analytical formula 

r,«(0, o)r+(-i,i)r+(-i, o)r ; ; (o. i) 

2w 



where w = 1^(0, o){r h "(o, i)[r+(-i, o) + r+ (-1, i)] + r+(-i, o)[r+ (-1, i) + r- fl (o, i)]} + r 6 -(o,i)[r+(-i,o) + 

r+ L (-l,0)]r+ L (-l,l). Here, we used R bL = R bR and denoted T+(-l,0) = 2r+ L (-l,0), 17(0,1) = 2^(0,1). The 
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auto-correlation function is derived from Eq. ( |14| 

q , n , r2 r t - fl (o,o)r+(-i,i)r+(-i,o)r b -(o > i) 

W 

In this case one see that the Schottky term S^£ h = 2SbL,bL(0)- The current through the top QD is 

T r- fl (o, o)r+ (-1, i){r- fl (o, i)r+(-i, o) + r~ (o, i)[r+ (-1, 0) + r+(-i, o)]} 

w 



(21) 



(22) 



The tunneling process between the states (0, 0) — > (— 1, 0), described by T ti? (0, 0), is relevant for I t L as well as SbL,bL- 
Therefore, the shot noise in the bottom QD (see Fig. 10) is proportional to the current in the top QD. 



C. Asymmetric resistances 



So far we have analyzed the system in the case when each QD is symmetrically coupled to the left and the right 
electrodes. In this section we will show how asymmetrical coupling to the leads can influence noises in the system. 
Let us analyze first the weak coupling case (Ci nt = 1 aF). The maps for the Fano factor F t L and the cross-correlation 
function StL,bL are presented in Fig. 11 The value of the Fano factor F t L is slightly smaller than in the symmetrical 



case [compare with Fig. [4^b)]. One can see that now the sub-Poissonian shot noise dominates in the region a. The 
factor F t L has non-monotonic behavior, which is well seen in the regions a, c and d. Its value is super-Poissonian in 
the area very close to the CB and close to the region e, next it drops below the unity with an increase of V t R and 
grows to the super-Poissonian values for higher voltages. The super-Poissonian shot noise occurs also in the regions e 
and g. We have found that, as in the symmetrical case, for an enhancement of F t L the auto-correlation processes are 
responsible, while the cross processes lead to a reduction of F t L below unity. The cross-correlation function S t L,bL is 
presented in Fig. 11 'b). In contrast to the symmetrical case [Fig. Qc)], when S t L,bL is always negative in the right- 
upper and the lefOower quarter of the map, now one see that it can be also positive in some areas. This dependence 
is non typical for cross-correlation in scattering of electrons and now it exhibits bunching. 

The results for the strong coupling case {Ci nt = 3 aF) are presented in Fig. 12 One see that behavior of F t L 



is very similar to that one obtained in the symmetrical case [compare Figs. |6jb) and (l2^a)]. The cross-correlation 
function S t L,bL in Fig- 12 'b) is non-monotonic, it can even change its sign with an increase of the bias voltage (see, 



for example, the region a and 6, where StL,bL < in the upper quarter and becomes positive for higher VtR). We 



made a cross-section for V t R 
to S t L,bL ~ see Fig. 



40 mV, through regions b and c, in order to analyze different contributions D 



n t n b ,n t n b 
tL,bL 
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For negative VbR the component D^®'® 1 is dominating (see the pink long dash curve). This 
contribution is also responsible for change of S t L,bL from positive to negative values with increasing VbR. When VbR > 



-11,01 



the amplitude of the component D t £\' L 10 and D t £^£ L is larger than the other components. Moreover, D 
describing correlations between tunneling events through the left top junction from the charge states ( 
forward tunneling processes in the bottom QD, becomes to play a dominating role. 



11,-10 
tL,bL ' 

1, 1) and 



V. SUMMARY AND FINAL REMARKS 



We performed studies of dynamical current correlations in the system of two capacitively coupled quantum dots 
(2QDs). In this 4-terminal system we first determined stability diagrams for states with various number of additional 
charge introduced either to the top or the bottom QD as a function of bias voltage V t R and VbR applied to the top and 
the bottom channel, respectively. The charge diagram is a background for studies current-current correlations, because 
we can separate allowed tunneling processes, which contribute to currents. We considered two different situations 
for weak and strong coupling between QDs. For weak coupling one can observe independent transport through one 
of QD in some range of bias voltages V t R and VbR (as one can expect). In contrast, for the strong coupling case, 
where both QDs are engaged in transport, the charging diagram is different. In the Coulomb blockade region the 
charge states ( — 1,1) or (1,-1) are metastable and when bias is larger than the threshold voltage they participate in 
transport [together with the state (0,0) and its neighbors]. For a very strong coupling a larger charge space with the 
states (— n, n) and (n, —n) (for n > 1) can be engaged. Transport through one of the channel induces large potential 
fluctuations in the second one [see, e.g., Sy ti y t and Sy hi y b in Fig.|5jb)]. This is a dynamical process and it is responsible 
for large enhancement of the Fano factor. We called it the dynamical Coulomb blockade effect) 29 We decomposed 

the correlation function Si a ^' a > into individual contributions D™^ n £ b /™, tnb of tunneling processes, which occur in the 
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charge space (n^n^). Dynamical contributions are relevant for the Fano factor. Their auto-correlated components 
are positive, whereas the cross-correlated components are negative. They compensate each other in part, and lead to 
the super- or the sub-Poissonian shot noise, depending on which contribution dominates (either auto-correlation or 
cross-correlation processes). 

A strong correlation can be seen even when one part of the device is unbiased, e.g., for V^r = and V t R > V th in 
the region b in Fig.^d). There occurs a pumping effect: charge transfer in the top QD leads to injection (or ejection) 
of charge in the bottom part of the device. Therefore, one can observe a large enhancement of the shot noise SbL,bL, m 
which the relevant contributions are from the injected I^ L and the ejected I^ L current t o the Schottky term Sjjj£ h . Our 
studies resemble shot noise in the presence of Coulomb drag in coupled quantum wires JSEU i n both the systems local 
charge and voltage fluctuations make electronic transport strongly correlated j 45 * 46 * which is seen in the current-current 
correlation functions. Of course our system is different, because we have assumed incoherent transport with electron 
thermalization at each quantum dot and electron transfer between the dots is prohibited, thus a momentum transfer 
and exchange processes are absent in our system - in contrast to coupled quantum wires, where these processes are 
relevant for the Coulomb drag effect. 

We showed also that the current cross-correlation function S t L,bL is in general negative, indicating on anti-bunching 
between charges transferred through the both channels. Let us stress again that here anti-bunching is a dynami- 
cal process caused by Coulombic repulsion between transferred electrons (strong potential anti-correl ations on th e 
both QDs), but it is not fermionic origin (as in recent studies for electrons scattered in nanostructures^ 12 * 13 * 14 * 1 ^ ) . 
The function StL,bL can have a plateau or can be non-monotonic as a function of bias voltage. The analysis of its 

components D^^ tnb showed that there are negative and positive contributions corresponding to forward and back- 
ward tunneling processes. Moreover, it was seen that due to retarded and advanced processes individual components 

^th'bL tnh = a ^ some finite bias voltages. All these various contributions give S t L,bL as a non-monotonic function, 
which can change its sign in some situations with asymmetric tunneling resistances. 

McClure et al. 1 performed an experiment on a 2QD system, which showed sign reversal of the current cross- 
correlation function S t ,b near a honeycomb vertex in the charge stability space. Applying gate voltages to the top 
and the bottom QD they could go from one stability region to another. Since the sign of S t ,b depends on directions of 
tunneling processes (whether electrons are injected to or ejected from QDs), sign reversal of St,b was observed at a state 
degeneracy line between the stability regions. The effect was also seen when one bias voltage was reversed. Notice that 
the experiment was performed close to border lines separating different Coulomb blockade regions and bias voltages 
were small. Our studies are more general and focus on dynamical aspects of current correlations, which occur for 
large bias voltages. It would be interesting to extend the experimental setupP^ outside the Coulomb blockade regime 
and verify our theoretical predictions, for example, to see competition of various contributions to the cross-correlation 
function S t L,bL- 

We presented also the volt age- volt age correlation functions Sv e ,v e ,, which are connected with charge fluctuations 
in our system. Since nanotechnology and shot noise measurement technique made recently great progress, we believe 
that such the correlation functions are not purely theoretical interest and they can be measured together with current 
correlation functions. 
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FIG. 1: (Color online) Schematic view of two capacitively coupled quantum dots. Each dot £ (£ = t,b — top, bottom) is 
connected with the electrodes (a = L, R - left, right). The tunnel junctions are characterized by the capacitances Ce a and the 
resistances Ri a . The strength of the dot-dot interaction is modeled by the capacitance dnt> The bias voltages VtR,VbR are 
applied to the right electrodes. 
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FIG. 2: (Color online) (a) The diagram of the available charge states and tunneling processes contributing to transport in 
the space of bias voltages (VtR, Vbii) for the coupling Ci n t = Ce a = 1 aF. Solid (green) lines denote threshold voltages V^, 
which set regions {a, . . . , n} with particular allowed tunneling processes {r| a (nt, rib)}, b) Diagrammatic representation of 
the tunneling processes {r| a (nt,nb)} between the charge states (nt,rib) allowed in transport for various regions {a, ...,n}. 
Tunneling processes between charge states are shown as (red) solid arrows for tunneling through the left junction and (blue) 
dash arrows denote tunneling through the right junction; horizontal arrows correspond to tunneling processes through the top 
QD, whereas vertical arrows - for tunneling through the bottom QD, respectively. For example, the diagram "6" presents 
allowed tunneling processes between states (0,0) and ( — 1,0) through the left and the right junction in the top QD. The solid 
(red) arrow between states (—1, 0) and (0, 0) denotes r^( — 1, 0), while the dash (blue) line between these states corresponds to 
F~ R (0,0). The other processes in the regions ...,n} can be described analogously. For larger voltages \V t R\ and (i.e., 
outside the regions enclosed by the green lines) additional new charge states appears, e.g., (2,0), (—2, 1), etc. 
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FIG. 3: (Color online) (a) The diagram of the available charge states and allowed tunneling processes contributing to transport 
in the space of bias voltages (VtR, VbR) for a strong coupling case d n t — 3 aF > CtL — CtR — CbL — CbR = 1 aF. (b) The 
diagrammatic representation of the tunneling processes \Xl a {n t , rib)} in the charge state space in various regions {a, . . . , n} 
(notation is the same as in Fig. |2|. 
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FIG. 4: (Color online) The current and the shot noise maps in the (VtR, VbR) space for the weak inter-dot coupling dnt = 1 
aF. (a) The current I t L [nA] flowing through the top QD. (b) The Fano factor F t L for the left junction in the top QD. (c) The 
cross-correlation function StL,bL(0) [10 _28 A 2 /Hz] for the currents in the top and the bottom QD. (d) The Fano factor F^l in 
the left bottom junction. The results were performed for R t L — RtR = 1 MQ, RbL — RbR = 50 MQ and T = K. The other 
parameters are the same as in Fig. [2] 
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FIG. 5: (Color online) (a) Plots of the Fano factor FtL (black solid curve) and all its components D t l t b £ t b contributing to 
the dynamical part of the shot noise StL,tL f° r VtR — 60 mV [the cross-section of the map in Fig. [4|b) through regions a, b, c 



and 

(-i 



The red dotted curve is for D. 



which corresponds for the auto-correlation function of the tunneling processes 



0) — > (0,0). Similarly, the blue long dash curve is for D~ L . The pink dash-dot-dot curve corresponds to D~ L ^ , 
which describes the cross-correlation of the tunneling events (—1,0) — > (0, 0) and ( — 1,1) — > (0,1). All the curves are normalized 
as the Fano factor (i.e., they are divided by 2eItL,)> (b) Plots of the voltage- voltage correlation function Sv t ,v t ( re d solid curve) 
and Sv b ,v b (blue short dash curve) describing potential and charge fluctuations at the top and the bottom QD, respectively. 
Note that the right hand axis corresponds to Sv b ,v b and potential fluctuations at the bottom QD are two order of magnitude 
larger than those in the top QD. 
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FIG. 6: (Color online) The current and the shot noise maps for the strong inter-dot coupling dnt = 3 aF. (a) The current 
ItL [nA] flowing through the top QD. (b) The Fano factor F t L for the left junction in the top QD. (c) The cross-correlation 
function StL,bL(0) [10 _28 A 2 /Hz] for the currents in the top and the bottom QD. (d) The Fano factor FbL in the left bottom 
junction. The results were performed for R t L = RtR = 1 MQ, RbL = RbR = 50 MQ and T = K. The other parameters are 
the same as in Fig. [3] 
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V tR [mV] 

FIG. 7: (Color online) (a) The Fano factor F t L (black solid line), F£ L (black short dash line) and its relevant components: 
D^ L 10 tL ~ 10 (red dotted curve), D^ 1 ^ 11 (blue long dash curve), and D^ ^ 11 (pink dash-dot-dot curve) plotted at VbR — mV. 
All the components are normalized as the Fano factor. The other parameters are the same as in Fig.p] (b) Voltage correlation 
functions Sv t ,v t (red solid curve) and Sv b ,v b (blue short dash curve) in the top and the bottom QD, respectively. 




FIG. 8: (Color online) The Fano factor FtL (black solid line) and its components: the Schottky term Ff L ch (black dash-dot 
curve) and the dynamical part F£ L (black short dash curve) for VbR = — 50 mV (a cross-section through the regions /, e, and 
d). The red dotted curve represents a component D~ L 10 tL ~ 10 , which is a relevant contribution to the dynamical part of the shot 
noise. The component D^l^ 10 is normalized as the Fano factor. In the region e the plot F>~ L 10 tL ~ 10 is very close to F£ L . The 
other parameters are the same as in Fig. [6] 
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FIG. 9: (Color online) (a) The cross-correlation function StL,bL (black solid curve) and the current IbL (black short dash curve) 
in the bottom QD plotted vs. VbR for the cross-section through regions a, b and c at VtR = 50 mV. The other parameters 
are the same as in Fig. 6 (b) The components of StL,bL'- the pink long dash curve corresponds to D^ ^® 1 , which describes 
correlations of the tunneling events ( — 1,0) — > (0,0) through the left top tunnel barrier and (0,1) — > (0,0) through the left 
bottom tunnel barrier; the red solid curve corresponds to D^ ^ 10 , which describes correlations between ( — 1,0) — > (0,0) and 
(— 1, 0) — > ( — 1, 1); the green dash-dot curve corresponds to D^^ 1 , which describes correlations between ( — 1, 1) — > (0, 1) and 
(0, 1) — > (0,0); and the blue short dash curve corresponds to D~^\^ 10 , which describes correlations between ( — 1, 1) — > (0, 1) 
and (-1,0) -> (-1,1). 
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FIG. 10: (Color online) The function SbL,bL(0) (solid black curve) plotted vs. the bias V t R for VbR = and d n t = 3. It is a 
cross-section through the region b of the map in Fig.|6jd). We also plotted its Schottky component Sb£ h (red dash-dot curve) 
and the frequency dependent part ^^^(O) (blue short dash curve). The plot shows that the shot noise in the bottom QD is 
induced by the current in the top QD I t L (pink long dash curve). 
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FIG. 12: (Color online) Maps (a) for the Fano factor F t L and (b) for the cross-correlation S t L,bL in the strong inter-dot case 
Cint = 3 aF and plotted for R t L = RtR = 1 MQ, RbL = 10 MQ, RbR = 100 MQ. The other parameters are the same as in 
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FIG. 13: (Color online) a) The cross-correlation function StL,bL (black solid curve) and the current IbL (black short dash curve) 
in the bottom QD plotted vs. VbR for the cross-section through regions c and b at VtR = 40 mV. The other parameters are the 
same as in Fig. 12 b) The components of StL,bL'- the pink long dash curve corresponds to D^®'® 1 , which describes correlations 



of the tunneling events ( — 1,0) — > (0,0) through the left top tunnel barrier and (0, 1) — > (0,0) through the left bottom tunnel 
barrier; the red solid curve corresponds to D~^ b ^ 10 and describes correlations between ( — 1,0) — > (0,0) and ( — 1,0) — > ( — 1, 1); 



the green dash-dot curve corresponds to D t 



and describes correlations between ( — 1,1) 



blue short dash curve corresponds to D t ^ 1 b , L 10 , which describes correlations between (—1,1) 



(0, 1) and (0, 1) - 
-> (0, 1) and (-1 



■> (0,0); and the 

o)-(-i,i). 



